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Assuming linearity of the perturbations at the time of decouphng, we reconstruct the primor- 
dial scalar potential from the temperature and polarization anisotropies in the cosmic microwave 
background radiation. In doing so we derive an optimal linear filter which, when operated on the 
spherical harmonic coefficients of the anisotropy maps, returns an estimate of the primordial scalar 
potential fluctuations in a spherical slice. The reconstruction is best in a thick shell around the 
decoupling epoch; the quality of the reconstruction depends on the redshift of the slice within this 
shell. With high quality maps of the temperature and polarization anisotropies it will be possible to 
obtain a reconstruction of potential fluctuation on scales between £ = 2 and £ ~ 300 at the redshift 
of decoupling, with some information about the three-dimensional shapes of the perturbations in a 
shell of width 250Mpc. 



I. INTRODUCTION 

Cosmic microwave background (CMB) radiation has 

S' ren us a great deal of insight into the early universe 
, for example the pattern of temperature fluctuation 
gives us the information about the perturbations in the 
plasma at a temperature T ~ 3000K. Using the concor- 
dance between big bang nucleosynthesis and CMB , we 
are able to probe directly the earlier universe T ~ 10^°K, 
but we rely on an understanding of cosmological theory 
to learn about the Universe at even higher energy scales. 
One of the most promising theories of the early universe 
is inflation 3], which predicts (i) spatial flatness of the 
observable universe today, (ii) nearly scale invariant, adi- 
abatic, primordial density perturbations, (iii) homogene- 
ity and isotropy on large angular scales of the observable 
universe, and (iv) Gaussianity. Predictions (i) and (ii) 
have been observationally tested while predictions 
(iii) and (iv) are yet to be conclusively tested. 

Any characterization of the Gaussianity (or otherwise) 
of the primordial perturbations will constrain models 
of inflation (e.g. [a, 0j 0)- Any non- Gaussianity pre- 
dicted in canonical inflation is very small 8, 9], but in- 
flation scenarios with large amounts of non-Gaussianity 
can be constructed [iflllllli3- ^^^^ bispectrum tests 
for primordial non-Gaussianity have not detected a sig- 
nal in the temperature fluctuations mapped by COBE 
[Hi and WMAP |[1J]. Other authors have found non- 
Gaussianity si gnat ures in the WMAP temperature data 
[Tsl ItR ITMll9j| . The main motivation of reconstruct- 
ing primordial potential perturbations is to study these 
non-Gaussianities. Reconstruction allows us to be more 
sensitive to the primordial perturbations, which is im- 
portant because current detections of non-Gaussianity 
do not specifically select for the primordial perturba- 
tions. An advantage of our approach is the ability to 



use the reconstructed maps to study higher order statis- 
tics of the primordial perturbations directly in the recon- 
structed potential. 

We generalize the results from to use both the 
temperature and polarization information. We will show 
that this results in a 3D reconstruction of a thick slice 
centered at Zdec- Using the polarization information re- 
moves the blind spots which are present at certain scales 
in the temperature-only analysis. We provide examples 
that show what will be possible with a future experiment 
that provide high-fidelity maps of the temperature and 
E-polarization anisotropies. 

We assume that the CMB two-point statistics (i.e. 
power spectrum) and other sources of cosmological in- 
formation will reveal the parameters to a high degree of 
accuracy such that they are input parameters in our re- 
construction. If desired our technique can be generalized 
to include residual uncertainties in cosmological param- 
eters by coupling it to Markov chains. 



II. RECONSTRUCTION OF PRIMORDIAL 
FLUCTUATION USING TEMPERATURE AND 
POLARIZATION 

The harmonic coefficients of the CMB anisotropy (aim) 
are related to primordial fiuctuation as: 



^frn = Z I k'dkr'dr[<S>,^{r) gf^{k)+S,Ur) 9xUk)]Mkr) 

(1) 

where ^emif) and Semif) are, respectively, the curva- 
ture and the isocurvature harmonic coefficients of the 
potential at a given comoving distance r = |r|; gxeif) 
is the radiation transfer function of either adiabatic or 
isocurvature perturbations; and X refers to either T or E. 
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Assuming that curvature perturbations $(r) domi- 
nate over isocurvature perturbations S'(f), we try to 
reconstruct the scalar potential $(r) from the observed 
temperature and polarization of the CMB. A linear 
filter Oe{r) which when operated on aim reconstructs 
the primordial field, can be obtained by minimizing the 
difference between the filtered field Oi{r)airn and the 
underlying field <^im{T) |23|. Since the perturbations 
were in the linear regime when the CMB decoupled, the 
choice of a linear filter is justified. Positive definiteness 



of 



7^ guarantees a minima of chi square at each 



and hence guarantees the existence of Of . 




d 



dOf{r) 



E Of{T)af^~^Ur)?)=Q (2) 



X=T,E 

which gives following form for the filters: 



Oj{r) 
Of{r) 



(jTE (jEE 



where 



(3) 



= («L«1™> = - / k'dk P^{k) gxt 9Yi{k) (4) 




Pfir) = {<S>irn{r)a}f,) ^-J k'dk P^{k) gxi{k) ,Mkr), 

(5) 

ji{kr) being the Bessel function of order £ and is 
the primordial power spectrum of scalar metric pertur- 
bations. In the presence of the instrumental noise, the 
expression for modifies. Assuming that the noise is 
Gaussian with variance af^, and spatially homogenious; 
we can write af^ 



c^QSu'Smm' then Cf-^ is replaced by (b)Of (r 



Cf^ + {(TqY- We have assumed that afj^ is zero. 
This defines an estimator for the potential 



X=T,E 



- (r) a 



(6) 



As evident the equations in (3) are coupled linear 
equations with variables Oj{r) and Of{r). 

If only one of the available CMB maps (either temper- 
ature T or polarization E) is used to reconstruct the 
potential, these equations decouple into ordinary linear 
equations in one variable Of{r) where X can either be 
T or E. In these special cases Eq. |21 modifies to 



^A\Of{r)af^-<^Ur)?) 



dOf{r) 

and gives the following form for the filter: 



(7) 



FIG. 1: Filters using T and E separately, at five different con- 
formal distances, r — c{to — Tdec), where To is the present con- 
formal time and Tdec is the conformal time at photon decou- 
pling epoch. We use a ACDM model with re-ionization and 

CTo = 14.2827Gpc and CTaec = 0.2836Gpc. (a)Oj(r) = 

I3f(r) 



We have first reconstructed the metric perturbations 
using the temperature and polarization information sep- 
arately as given by Eq. |S1 Then we solved Eq. El to 
combine the information from temperature T and E po- 
larization and show that reconstruction is much better. 

Using both T and E, the expected variance of the re- 
constructed potential is given by 



where aj^ and af^ are the noise variances in the tem- 
perature map and in the polarization map respectively. 
We define 



Ot{r) 



(8) 



(10) 
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and use its expectation value: 

(Pf ) = i0jir)rcr+i0f{r)rCf^+20j{r)0f{r)Cr, 

(11) 

for the purpose of comparing our reconstruction to the 
theoretical power spectrum. As an estimator of the pri- 
mordial power spectrum this biased because our filters 
are optimized for estimating the potential ^imif) and 
not the power spectrum. This is the standard for opti- 
mal, or Wiener, filters, since they remove contamination 
aggressively. One can compute the bias to correct for this 
effect and bound the allowed range for the true power 

spectrum as ij^f'^ i ^li where 

= V2{p/(r) + {Oj{r)fCr + {Of{r}fCr 
-2(0f (r)/3j(r) +Of (r)/3f (r) - Of(r)Of (0/3^(0) }{I2) 

Here, Pf{r) is the theory power spectrum of the potential 
on a spherical shell of radius r. 



not minimize the error bars on this power spectrum. The 
fact that the error bars include the true spectrum show 
that the reconstruction of the perturbation itself is unbi- 
ased. At small angular scales t ~ 300 subhorizon physics 
wash out the information about the primordial fluctua- 
tions and the reconstruction fails. 

The combination of temperature and E polarization 
data enables the reconstruction up to ^ ~ 300, without 
blind spots in t. This is a marked improvement compared 
to the performance of the operator that uses either tem- 
perature or polarization alone which are shown in Figures 
3f and 4f. 

In Figures 3 and 4 the error bars increase due to cosmic 
variance on large angular scales. In Figure 3 we see that 
re-ionization improves the ability to reconstruct the pri- 
mordial potential at large scales (except in a trough near 
£ ^ 20). The polarization-only reconstruction reveals 
that the enlarged E polarization signal due to early re- 
ionization aids in the reconstruction at very low I. Physi- 
cally, the rescattering of photons at low redshift imprints 
information on the large scale modes of the primordial 
perturbation . 



III. RESULTS 

Our examples are computed within a ACDM model 
with and without re-ionization (Oj, = 0.26, fib = 

0. 04, r^A = 0.7, h = 0.7, scalar spectral index Us = 0.93, 
with running — —0.031, re-ionization optical depth 
T = 0.17). 

Figure 1 shows plots of the operators Of{r) (for X 
as T and E respectively) for different r as a function of 

1. While the reconstruction of the primordial potential 
will be good on large scales for both T and E (agreeing 
with the Sachs- Wolfe effect ^l]) , it will be bad on small 
scales as Of{r) oscillates about for both T and E. But 
as a consequence of the physics at decoupling, the trans- 
fer functions Qt and gEi{k) are out of phase. This 
results in Oj{r) and Of{r) never both being zero at a 
particular scale (see Figures 3f and 4f). This observation 
supports our analysis leading to the coupled Eq. |21 and 
we are able to reconstruct the potential on much broader 
range of scales if we use both temperature and polariza- 
tion information from CMB. 

In Figure 2 we show reconstructed potential maps at 
various slices of conformal time. The quality of recon- 
struction is sensitive to the reconstruction depth and is 
best in the thin slice centered on the peak of the dif- 
ferential optical depth. By stacking these slices a three- 
dimensional thick shell of primordial potential is recon- 
structed. 

Figures 3 and 4 show the expected power spectra of the 
reconstructed potential maps for models with and with- 
out re-ionization. The power spectra are bounded from 
above by the theoretical power spectrum computed from 
our cosmological parameters Ug and A. The amount of 
reconstructed power is indicative of the quality of the re- 
construction. As discussed in section II, our filter does 



IV. CONCLUSION AND FUTURE WORK 

We have demonstrated the application of linear the- 
ory to construct a least-square estimator of the primor- 
dial potential of adiabatic density fluctuations generated 
during the epoch of inflation. Starting with high-quality 
maps of the CMB anisotropics in temperature and E po- 
larization and using a fiducial set of cosmological param- 
eters we can reconstruct the primordial potential in a 
spherical shell of thickness ^ 250Mpc near the peak of 
the differential optical depth. Reconstruction is reason- 
ably good up to £ '--^ 300 and the stronger E polarization 
signal due to re-ionization aids in reconstructing the po- 
tential at very low £. This effect is consistent with and a 
generalization of the results by Skordis and Silk _22 j who 
reconstruct the CMB quadrupole with accuracy better 
than cosmic variance. 

To exploit the tomographic techniques we have devel- 
oped, we require maps of T and E and high signal-to- 
noise ratios. For this first exploration we have neglected 
the effect of foregrounds, which will reduce the effective 
area of the T and E maps that we will be able to use in 
the reconstruction. One of the advantages of CMB anal- 
ysis techniques based on Gibbs sampling 01 is that 
methods such as this one can be applied essentially with- 
out modification while the messy details of dealing with 
systematics like foregrounds and instrumental contami- 
nants are taken into account by the Bayesian analysis. 

One might wonder whether there are ways to improve 
on our techniques. In principle this could be done in 
two distinct ways. First, we might wish to broaden the 
range in redshift over which we can reconstruct the pri- 
mordial perturbations. A second way to improve on our 
techniques would be to enhance the quality of reconstruc- 
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FIG. 2: Reconstructed potential maps in slices centered around the sphere of last scattering r^ec = c{to — Tdec)- The maps 
are ordered in decreasing distance, left-to-right and top-to-bottom. The color scale is in dimensionless units of fractional 
perturbation. Each map is 25 degrees on the side. 



tion. 

On the first issue, scattering of cosmic microwave back- 
ground photons through processes other than Thompson 
scattering, such as Rayleigh scattering or resonant scat- 
tering makes the location of the last scattering surface 
frequency dependent[2^ In principle, this type of 

effect would allow a more detailed tomography of the in- 
flationary perturbations over a broader range in redshift. 
However, in the models currently favored by cosmological 
observations Rayleigh scattering only varies the redshift 
of the peak in the visibility function by about 1% [25| . 
This is much less than the thickness of the last scattering 
surface, and hence does not broaden the redshift range 
where our reconstruction is effective. A detailed study of 
the promise of using resonant scattering for this purpose 
remains to be done. 



On the second issue, improving the quality of the re- 
construction, it is possible that additional linear combi- 
nations of primordial perturbation modes could be con- 
strained by high signal-to-noise CMB maps over a range 
of frequencies. 

Ultimately we envision tests of non-Gaussianity ap- 
plied to the reconstructed primordial potential. Our fil- 
ters will select those combinations of the spherical har- 
monic modes that correlate to the primordial potential. 
This reduces the probability that residual foregrounds 
result in a non-Gaussian signature. Further, the ability 
to reconstruct the potential in a thick shell allows the 
application of three dimensional shape statistics to po- 
tential peaks and troughs. We anticipate the search for 
non-Gaussianities in the primordial perturbations once 
high resolution and high sensitivity maps of the E polar- 
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FIG. 3: Power spectrum of the primordial fluctuation for the ACDM model with re-ionization {Qc = 0.26, Qt = 0.04, Qa = 
0.7, h = 0.7, T = 0.17, Us — 0.93, — —0.031). For all plots the solid line shows the theoretical primordial power spectrum. 
In (a)-(e) we show the reconstructed power spectrum using the coupled estimator given by Eq. (|SJ (dashed), bracketed by error 
bounds (dot-dashed). In (f) we show the reconstructions due to temperature (dashed) and polarization (dot-dashed) alone, 
using the estimator given by Eq. ||5J. The theory and the reconstructed power spectrum are shown in arbitrary units. 



ization signal become available. To our knowledge the 
techniques we describe in this paper represent the most 
direct way to date to probe the statistical properties of 
the primordial perturbations created in the inflation era. 
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FIG. 4: Power spectrum of the primordial fluctuation for ACDM model without re-ionization (Sic = 0.26, Q,h = 0.04, S7a = 
0.7, h — 0.7, r = OjTis = 0.93, — —0.031). For all plots the solid line shows the theoretical primordial power spectrum. In 
(a)-(e) we show the reconstructed power spectrum using the coupled estimator given by Eq. © (dashed), bracketed by error 
bounds (dot-dashed). In (f) we show the reconstructions due to temperature (dashed) and polarization (dot-dashed) alone, 
using the estimator given by Eq. ||5J. The theory and the reconstructed power spectrum are shown in arbitrary units. 
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